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I. 



HWRODUCTION 



With much renewed interest in transonic flow, the limitations of 
present methods of analysis become evident. The hodograph transformation 
can be useful only for flows past a small class of airfoils; numerical 
methods are powerful but lack some generality. It is, therefore, worth- 
while to examine anew the possibility and suitability of exact solutions 
to the two-dimensional, small disturbance potential equation of transonic 
flow. This equation retains one non-linear term which is essential for 
non-divergent solutions at Mach one, but the non-linear partial differen- 
tial equation has proven to be difficult to solve. This inherent difficulty, 
coupled with the presence of shocks in the flow which cause boundary layer 
separation, has resulted in the creation of many approximate methods of 
solution^ which are employed in the design of transonic airfoils and the 
like . 

Some exact solutions are available for the transonic equation which 

Q 

are not obtained through the hodograph method. Mayer obtained a solution 
for the flow in a deLaval nozzle in the form of a series expansion; 

Goertler and Guderley^ report on an exact solution to the "parallel sonic 
jet." Exact solutions are limited either by boundary conditions as in 
the hodograph method or by the inherent approximations of the transonic 
equation. These solutions, however, can be very useful in checking 
results from approximations to the equation and do represent a practical 
class of geometries in their own right. 

In this paper an exact solution is presented which satisfies the full 
form of the transonic flow equation and involves 5 parameters . This solu- 
tion has been obtained through a separation of variables which is appro- 
priate to a re-arranged form of the original equation. These resvilts can 

be studied in the phase plane where significant features of the transonic 
equation are evident. The inipllcit form of the solution can be given an 
explicit asymptotic form which is both simple and useful. The separated 
solutions can be manipulated to generate transonic, shockless flows. 
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II. SEPARATION OF VARIABLES 



3 

The transonic equation in two dimensions is 
where 

u . V. 

'P^ "" tT 'Pv " iT 

X Uqq y uqo 

u and V are the perturbation velocities corresponding to the x and y axis 
respectively, and is the Mach nuniber of the unperturbed flow at Uqo 
moving in the positive x direction. The coordinates are dimensionless. 

Equation 1 can be separated with 



cp(x,y) = X(x)Y(y) + (l-lV^)x 

where 

^(x,y) = m| (1+y>P 



( 2 ) 

(3) 



the resulting form is 



X*X” 

X 




- X 



(X>0) 



where \ is a separation constant the sign of which is chosen so that Y 
remains finite as y goes to infinity. In the subsonic wavy wall solution^ 
it is shown that the separation constant is related to the "wavelength" of 
the sinusoidal wall. We arrive then at two ordinary, nonlinear equations 



o 

II 

+ 


(4a) 


+ XY^ = 0 


(4b) 



To obtain the X(x) solution, let P(x) = X'; integrating once we obtain 



P(x) 



|xX^ 



(5) 



2 



where a is a constant of integration. Now, integrating again we get 



X 






X 




dZ 




+a 



(6) 



which is an implicit form for the solution. The Y(y) solution follows in 
a similar fashion; letting Q(y) = Y' we arrive at 




2 3 

^Y^+P 



( 7 ) 



where P is an integration constant. Upon integrating again the result is 



y = Yo 




(8) 



which represents the inplicit solution for Y(y). 



We shall study next the separated solutions in the phase plane 
before further considerations of the above integrals are given. 
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III. PHASE PLANE 



Since the independent variable does not appear in Eqs. 4a and 4b, 

4 5 

then the analysis in the phase plane is appropriate Using the previous 
definitions for P and Q we can rewrite Eqs. 5 and 7 as 





ip3 . = i ./3 


(9) 


and 


= + p/2 


(10) 



Provided that neither a nor p are zero, it is useful at this point 
to divide Eqs. 9 and 10 by a and p respectively. What results is an 
equation of the form 





a^ + b^ = + 1 


(11) 




P /2 


(12) 


where 




and 


b = (| X or ^2 


(13) 



It is a simple matter to plot Eq. 11 on the phase plane and this is 
shown in Fig. 1 . In this figure, we have also shown the corresponding 
curves for a circle ( a degenerate ellipse), and for a hyperbola. Since 
both the ellipse and the hyperbola represent solutions to linear equations, 
the non-linear character of the present solution is evident in Fig. 1 . 

In order to further explore the features of the phase plane, let us 
return to the velocity potential; from Eq. 2, 

cp^ = X*Y + (1 -mJ) (14) 

and S = XY’ (15) 

*y 



where 



u 



UoomJ(i-»y) 



and cp = 



U, 



fco^£(l-^v) 
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FIGURE 1. HiASE PLANE REPRESENTATION 
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Rearranging Eqs. l4 and 15 



X'Y = cp 






(16) 



xr = cp 



(17) 



Clearly, a phase plane representation of X' (or P) vs X corresponds to 

plotting [cp^ - (l-]^)] vs cpy with a constant y coordinate. What is equally 

significant is that the same figure results if we plot Y vs Y' (or Q) along 

a constant x coordinate. The above inference can be seen from Eqs. 11-13. 

Figure 2 shows the phase plane in cp and cp ; Z signifies either Y or X' 

^ y 

and W signifies either Y* or X depending on which coordinate is being held 
constant. The character of Fig. 2 is useful because our inquiries as to 
equilibrium and stability points in the phase plane will have direct 
physical significance. We shall also find that trajectories will be of 
some interest. 

It can be seen in Fig. 2 that, except for = 1, the curves do not 
go through the origin. Since we know that there must exist a region for 
which no velocity perturbations exist (i.e. cp = c^ = O), it is clear that 
intercepts in the ordinate are proportional to 1 -m£ . That is, the upper 
curve must correspond to an initially supersonic flow (l-M^ o) and the 
lower curve to an initially subsonic flow; the axis is the no pertiirbation 
point. Although shocks are not necessarily normal but must satisfy the tran- 
sonic shock polar, the above argument is consistent with normal shocks near 
Mach one where 

m |^-1 = 1 -M^ (18) 

holds between supersonic (subscript 1) and subsonic Mach numbers. Clearly, 
without pertiirbations , a flow may shock and this corresponds to a trajectory 
along the ordinate between the upper and lower curves. 

At this point we are ready to interpret 0/ and p since their sign 
governs the choice of curve in the phase plane . It follows that 



and 



+a=C^ (M^-1) 



+ 3 (M^-1) 



(19) 

( 20 ) 



6 



(I-M^ 





FIGURE 2. FERTUEBATION VELOCITIES IN PHASE PLANE 
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where and are constants. The form of Eqs. 19 and 20 can he shown 
to he consistent with the von Karman's transonic similarity rule as 

O 

given by Spreiter"^ (where for a fixed value of x and y, X'Y must be 
a function of see Section V). 

Using phase -plane arguments* it follows that the origin in Fig. 2 is an 
equilibrium point for This means that, at least for two dimensions, 

such a flow cannot develop any perturbations since for P and X equal to , 
zero there is no way the curves may leave the origin. 

It is instructive to look at the trajectories in the P-X plane and 
Q-Y plane separately. The latter are shown in Fig. 3- The arrows indi- 
cate the direction of increasing y. The origin is an equilibrium point 
(Q^=0 and Y’=0) which iniplies that the only "stable" trajectory is that 
for Moo=l. Q = 0 corresponds to an ‘asymptote; because y->co as Q-»0 from above 
and y-*-oo as Qr*0 from below; however, Q always decreases for increasing y. 

The P-X trajectories are more coniplicated as shown in Fig. 4 . The 
trajectories for P>0 and R:0 are opposite iicplying that the P=0 intercepts 
are " singular "point s ; note, however, that these are points at x=+oo. The 
subsonic side (lower curve) can be traversed without difficulty but the 
supersonic side shows problems in going from region I (P<0) to region II 
(P>0) and viceversa. This difficulty will be easily understood in Section 
IV where curves are shown in the coordinate plane. 

The P-X plane shows trajectories which may elucidate transonic shock 
behavior. In Fig. 2 the perturbation velocities are indicated for the motion 
represented by Fig. 4 (i.e. with x increasing) and for Q > 0. A supersonic 
flow would start at the ordinate intercept of the upper curve (cp^ = cp =0) 
and move to the right whereas an originally subsonic flow would begin at the 
lower intercept and move to the left. A variety of flows are possible depend- 
ing on the choice of the free parameters involved; however, the choice is re- 
stricted to region II for initially supersonic flow until we introduce transonic 



*The equation in X(x) must be modified as follows 

let R * I (X')^ 
then X' =V^ and R' = - \X 
8 




FIGUEE 3. Q-Y TRAJECTORIES, y INCREASING 
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X 




FIGURE k. P-X TRAJECTORIES, x INCREASING 
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shocks. The discussion of a subsonic flow which goes through a local 

supersonic zone within the airfoil appears to be more complicated and will 

not be attempted here. The important question as to whether or not shocks 

can be avoided in transonic flow remains, to be shown with the phase plane analysis. 

IV. IMPLICIT SOLUTIONS 

In order to present further discussion of the implicit solutions, 
we note that for appropriately large values of X and Y the following 
simplifications take place: 

P(x) - |xx^« - - 

and Q(y) (|x)VV/= 

where only positive values of Q have been chosen. 

This is, of course, evident from Figs. 1-4 where it can be seen 
that the subsonic and supersonic branches approach a common asymptote so 
that at sufficiently high values there is no distinction between the sub- 
sonic and supersonic regimes within transonic flow. The approximation, 
furthermore, represents the exact solution for Mo,=l. 

Under the above conditions, the integrals represented by Eqs. 6 and 
8 become 

X = - ^ x^ (21) 

and ^ " ■ f ^^"^ 0 ^”^ 

(Note that Y = -“aty = y^, whereas X=^0 at x=0) 

Therefore, from Eq. 2, 

?U,y) = I — — 2 + (23) 

3 (y-y„)" 

where, \, the separation constant related to a geometry input has dropped 
out. The limit represented by Eqs. 21 and 22 is suggestive of supersonic 
flow at regions aWay from the body surface. Consider the perturbation velo- 
cities , 2 

5 — o + (s'*) 



xx|)V^ 

Y<-(|)^/= 



11 



and 



(25) 



= -3 



(y-yj- 



we 



If we "force^* the characterics for supersonic flow on these results, 
find that, using x = + (M^-1)^ (y-y ), Eqs. 24 and 25 become 



'Jp = 0 (26) 

X 

^y " ^ 3 (27) 



It is interesting to compare these results with a Prandtl-Meyer 
expansion for values close to the sonic line. The first non-zero term 

/ N 6 

of the angle V (which changes M=1 to H;o>l) is , using a series expansion 
for the tan 



_ 2 (£-1) 

" 3 (Y+l) •" 




(28) 



but from Eqs. 26, 27, and 3 



V ^ tanv 



V 

Uoo-Hi 



2 (m£-1)^^^ 

3 (y+1)H»^ 



(29) 



The comparison is good since Moo niust be close to 1. If Moo is exactly one, 
then we know that v =0; the sonic solution is simply 



~ 1 
) = 5 



{y-y^f 



(30) 



It is evident that if x is zero V , 'P , and cp will all be zero which 

x’ xx’ yy 

corresponds to the trivial solution to the transonic equation. This is 
the equilibrium point alluded to in the previous Section. 

2 

The parallel sonic Jet solution of Goertler and Guderley is con^rised 
of Eqs. 8 and 21, but note that sign of \ as chosen by these authors is 
opposite that used in this work. The result is to make their Y(y) go to 
infinity at some finite y, and this appears to be inconsistent with the 
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character of unbounded flow. As will be shown at the end of this Section, 
our Y(y) has an upper limit as y goes to infinity. It is possible, 
however, to obtain Eq. 30 from the Goertler and Guderley solution. 

Before proceeding to depict the entire range of the integrals given 
by Eqs. 6 and 8, it is worthwhile to examine their behavior near X=0 and 
Y=0. We have, using Eqs. 19 and 20 as well. 



which says that the slope of X near the origin is constant, positive for 
supersonic flow and negative for subsonic flow. In the y direction we 
have 



which implies that only the supersonic flow solution crosses the abscissa 
and with a positive, constant slope. Because of the third power on Y, 
the straight line portion of it will be conspicuous. 

If we now write Eqs. 6 and 8 for Moo^l as 




(31) 




(32) 




(33) 





and define 



X 3 X O' 






(35) 




(36) 



13 



then we can plot X vs 3c' and Y vs (y-y^). This is shown in Figs. 5 and 6. 
These results were obtained by numerical integration of the equations. 

In Figure 5> the axis was chosen as the point of symmetry since the curves 
go to infinity. It is clear that the supersonic solutions of region I 
cannot pass into region II (except at +<„) as ircplied by the phase plane 
plot, Fig. 4. The subsonic solution has no discontinuities. For values 
of X greater than + 10, the form of Eq. 21 well represents the curves 
within three significant figures. In Fig. 6 we have plotted only the 
curves for the curves for would correspond to a reflection of the 

shown curves about the ordinate. Since Q, always decreases with increasing 
y, crossing the abscissa in Fig. 3 requires a trajectory going through 
infinity. It is clear in Fig. 6 that only the supersonic solution crosses 
the y-axis (at y-y = + 2.8) and that it is asymptotic to Y=1.0; similarly 
the subsonic solution is asymptotic to Y=-1.0. For values of Y less than 
-5 the curves are well represented by Eq. 22. With knowledge of X(x) and 
Y(y), P and Q, can be found from Eqs. 5 and 7 respectively. 

V. INTRIWSIC FLCWS 

The solutions which arise from separation of variables in the transonic 
equation generate three flows which we shall call intrinsic. The reasoning 
will be evident in the arguments shown below. Eqs. 33 and 3*+ may be written 
as (M^?a) 







Now at the body surface y.^^ = f (x) or 

\ = g (X) (39) 



Let us define the intrinsic flows as ones for which 



Y, ^ = -X ±1 

b 



(40) 



It is not difficult to show that as a result of Eq. 40, taking derivatives. 



dY, 






2 

3 



d^ 

r ^2+ , 2/3 

L-X - Ij 



(41) 
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X 




FIGURE 5. NUMERICAL INTEGRATION OF EQUATION 33 

X5 



Y M3,>i 




FIGURE 6 . NUMERICAL INTEGRATION OF EQUATION 3^ , 
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( 42 ) 



Introducting Eqs. 37 and 38? the above becomes 



dx 



2 

3 






Since, according to Figs 4 and 5> there are three distinct curves (one 
subsonic and two supersonic) we would anticipate that Eq. 42 represents 
three flews. Before going to a numerical integration, let us see what 
Eq, 42 yields for both the X > 1 and X < 1 limits. We have 



and 






constant 


- 6/s. 


for X> 1 


(43) 


3 


X < 1 


supersonic 


(44a) 


L 2 w 

(+3 X, 


X < 1 


subsonic 


(44b) 


coefficient in 


two-dimensional flow at the 


surface 


3 








iPa-UoD^ 


= - 


2 (^.y^) 


(45) 



be found as follows 

C ^ 
P 



But, according to Eqs 3 and l4 _ 

2[x'y^ + (1-mJ)] 

% (1 + Y) 



(46) 



It can further be shown that 



X'Y. = [ llf/3 (c^/3) (» 



1/3 



2X' 



where ( q ? 0)^'^3 = (c.C_)^'^3 (m ^- l )^'^^ 

± d 00 



(4?) 

(48) 



So that C becomes 
P 



C = - 



Mfl, ( 1 ^) 






) ^(49) 



Eq. 49 shows the dependance of C on (l-M_^) explicitly. We next define 

P 

. __ _ ,/3 

S 2/3O1C2U/3 (M 2 .i) 2/3 



/V/ 



( 50 ) 



where is the part of the pressure coefficient which varies with the 
coordinate of the surface. 
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Of the three possible flows, only two turn out to he reasonable, namely, 
the one corresponding to M^< 1 and that for M^> 1, region I. These flows 



00 



are shown in Figs 7 and 8 where y^ and C are plotted as a function of 3c. 

The initially subsonic flow represented in Fig. 7 corresponds to a concave 

shape connected to a convex one through a strai^t section; there is no 

stagnation point. The curvature of C reflects the concave -convex features 

of the surface. The maximum value of C is minus one at x = 0. We note 

P 

that, since Moo< the conventional C will always be negative unless X 
itself is not as previously defined. The initially supersonic flow shown 
in Fig. 8 is to be interpreted as flow in the minus x and minus y directions 
(i.e., to the left and dcwn). Here the shape is all concave and we expect 
few differences between this flow and the previous one. C , for example, 
goes to zero at x = 0 where one would anticipate noticeable differences as 
evidenced in Fig. 5. These results will, of course, break down at some y^ 
where the small perturbation assumptions begin to be exceeded. These two 
flows emerge in a straightforward way from the form of the separated 
solutions and should become part of the important class of transonic flows. 
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FIGURE 7. INTRINSIC FLOW, SUBSONIC INITIALLY 
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FIGURE 8. INTRINSIC FLOW, SUPERSONIC INITIALLY REGION I 



VI. COIIDLUSIOI^ 



The ultility of the exact solution shown in this paper is associated 
with cases where the boundary conditions are compatible with the separa- 
tion of variables scheme. This is particularly significant on the body 
surface (the so-called tangency condition); a family of flows which is 
intrinsic to the form of tbs solutions has been discussed. At. 
infinity, the perturbation must be finite or zero and the fact that Y(y) 
is bounded as y-»oo means that this boundary condition may be handled as 
well. Of course, there is a class of internal flow problems for which 
our solution may also apply. 

Since the phase plot may be interpreted in terms of the perturbation 

velocities and cp , we have before us a fundamental tool for studying 
X y 

the general properties of the transonic equation. Paths of constant x or 
y coordinate are easy to see in practise and one can find trajectories 
for which the pertiirbation velocities will remain constant. The allowed 
trajectories with increasing x promise to be particularly useful in un- 
covering the behavior of shocks in transonic flow. Moreover, the conclu- 
sions from the phase plane that govern sonic flow point out some rather 
fundamental behavior of two-dimensional flow, namely, that sonic flow is 
an equilibrium point. 

The asymptotic solution to the velocity potential has been identified 
with sonic flow and with a very weak Prandtl-Meyer expansion. Certainly, 
the siirplicity of the mathematical form of this result makes it worthy 
of f further study. 

The numerical integration results given have aided in explaining 
some of the features of the given solutions. The characteristics of the 
flows presented will hopefully have an iit5)act on practical designs , 
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